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ABSTRACT 

We present RXTE monitoring of the eclipsing X-ray binary Hercules X-l conducted over 
the short-high state of July 1998. This was one of the last major short-high states before the 
source entered an anomalous low-state of activity. A comparison with previous epochs finds no 
evidence for special behavior during these observations. We determine orbital and pulsar spin 
periods to facilitate measurements of P spin and P arb during the subsequent anomalous low state 
and during the next epoch of high-state activity. Spectrally, the decay of the short- high state and 
concurrent pre-eclipse dips are consistent with obscuration of a central X-ray source by a cloud 
of non-uniform column density. The standard model of a warped accretion disk of finite vertical 
scale height fits the characteristics of this absorber well. Pre-eclipse dips have durations a factor 
of a few longer than the characteristic durations of dips during main-high states. Pulse profile 
structure increases in complexity towards the tail of the short-high state suggesting changes in 
accretion curtain geometry. 

Subject headings: accretion: accretion disks - binaries: close - binaries: eclipsing - stars: individual 
(Hercules X-l) - stars: neutron - X-rays: stars 



1. INTRODUCTION 

Her X-l (Tananbaum et al. 1972) is an eclips- 
ing X-ray binary containing a pulsar of 1.4 M 
and an A7 stellar companion of 2.2 M Q (Middlcd- 
itch & Nelson 1976; Reynolds ct al. 1997). The 
system displays behavior on four separate period- 
icities - the pulsar spin period (1.24-s), the bi- 
nary orbit (1.7-d), a super period of 35-d which 
results from a retrograde-precessing, warped ac- 
cretion disk, and a beat between the precessional 
and orbital periods of 1.62-d. The warp engine is 
poorly understood but is likely to be a combina- 
tion of both radiation- and tidally-driven preces- 
sion (Papaloizou & Terquem 1995; Pringle 1996). 

In X-rays the 35-d super-cycle results in two 
phases of strong X-ray activity per cycle (Giac- 
coni et al. 1973; Scott & Leahy 1999). The main- 
high has a rapid turn-on over ~90-m and decays 
over ~10-d. A low-state where flux remains at 3-5 
percent of the high-state level follows for the next 
<~10-d, succeeded by a second short-high state, 
lasting ~5-d, with flux peaking at ~30 percent 
of the main-high state maximum. A further low- 
state completes the cycle and extends over the 
next ~10-d. UV, optical and infrared wavelengths 
are dominated by X-ray emission reprocessed over 
the inner face of the companion star. There is no 
high-low state switching in these wavebands be- 
cause the pulsar is always visible to a large frac- 
tion of the companion star's surface, however the 
accretion disk casts a shadow over the stellar pho- 
tosphere, whose terminator migrates on the pre- 
cession timescale. Consequently, the super-cycle 
is also observed at UV, optical and infra-red ener- 
gies (Gerend & Boynton 1976). 

The 35-d clock has remained quasi-coherent 
since it's discovery, although there have been two 
recorded occasions when the clock has missed 
several consecutive turn-ons (Parmar et al. 1985; 
1999), and one occasion when the X-ray flux dur- 
ing main-high states was significantly reduced (Vr- 
tilek et al. 1994). The cause of these anomalous 
low-states is not clear, but they probably result 
from changes in the state of the accretion disk, ei- 
ther through an increase in vertical scale height, or 
modifications to the disk warp. However the rela- 
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tively constant level of UV and optical flux during 
these states suggests that accretion does not turn 
off altogether. 

During spring 1999, Her X-l entered its latest 
anomalous low state and remained there until the 
fall of 2000. Fig. 1 presents the long-term light 
curve of Her X-l, sampled once per day over the 
energy band 2-10 keV, as obtained by the RXTE 
All-Sky Monitor (ASM) on board RXTE (courtesy 
of the ASM team; Levine et al. 1996). After MJD 
51200 there are two weak main- high states before 
the onset of the anomalous low-state. The data 
presented in this paper were taken at MJD 51000 
which appears to be a normal short-high state. 
However all subsequent short-high states appear 
either to be absent or weak within the sensitivity 
limits of the ASM. This suggests that a precursor 
to the anomalous low state may be a weakening of 
short-high states beginning some 200 days before 
the main-high states vanish. This would have con- 
sequences for the timescale of mass transfer vari- 
ability in the system, and the dynamical reaction 
time of the disk warp to mass transfer rate onto 
the neutron star. 

The same short-high state was simultane- 
ously observed with BeppoSAX. Oostcrbrock 
ct al. (2000) present data obtained with the Low- 
Energy Concentrator Spectrometer (LECS; 0.1-10 
keV) and the Medium-Energy Concentrator Spec- 
trometer (MECS; 1.8-10 keV) experiments over 
four orbital cycles. They discern long-duration, 
energy-dependent dips in the light curve, model 
spectra as composites of powerlaw, blackbody, and 
Fe line components, provide energy-resolved pulse 
profiles and fit pulse-phased spectral models to the 
data. The 100 eV blackbody component displays 
pulses which are out of phase with the powerlaw 
component. It is believed to originate from the 
inner regions of the accretion disk. 

Oosterbroek et al. (2000) provide spectral 
models and spin-phased light curves for the least- 
absorbed intervals of each orbit, but do not follow 
the spectral or pulse profile evolution through 
dipping events. RXTE's larger collecting area 
compared with BeppoSAX allows us to present 
pulse profiles and spectral fits through the dipping 
events, and over an energy band harder than the 
BeppoSAX LECS and MECS instruments. In the 
following paper we also search for any outstanding 
properties of this short-high state that might sig- 
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Fig. 1.— The RXTE ASM light curve of Her X-l. The label "PCA" indicates the time and combined 
duration of the current short-high state visits. 



nal the onset of any future anomalous low-states, 
and measure the pulsar orbital and spin periods 
to help enable the correct measurement of orbital 
and rotational changes during the low state. These 
data form part of a simultaneous multi- wavelength 
campaign of Her X-l(Vrtilek et al. 2001; Boro- 
son et al. 2000; Boroson et al. 2001; this pa- 
per and Still et al. 2001), consisting of pointings 
from RXTE, the Extreme Ultra- Violet Explorer 
(EUVE), the Hubble Space Telescope (HST), the 
4.2-m William Herschel Telescope on La Palma 
and the 3.5-m Calar Alto telescope, Spain. 

2. OBSERVATIONS 

RXTE pointed at Her X-l intermittently dur- 
ing 1998 July 9-13 (MJD 51003.8-51007.9) accu- 
mulating 70-ksec of time-tagged events. Assuming 
the 35-d epoch is determined by (M. Kunz, private 
communication) : 

T35 = MJD 50041.0 + 34.85 E35 (1) 

where T 35 corresponds to X-ray turn-on of the 
main-high state and E 35 is the cycle number, the 



precession phases sampled are 35 = 0.63-0.74. 
This range corresponds approximately to the peak 
of a short-high state to the end of the short-high 
state, respectively (Scott & Leahy 1999). 

Event reconstruction was performed using stan- 
dard algorithms within ftools v4.2. We anal- 
yse data from the Proportional Counter Array 
(PCA; Zhang et al. 1993) which consists of five 
identical Xe Proportional Counting Units (PCUs) 
with a combined effective area of 6500 cm 2 . As 
well as the standard data formats, we employ two 
event analysers in the GoodXenon event mode 
with 2-s readout. Timestamps are resolved to 1- 
fis while employing the full 256 energy channels 
(1.8-101.0 keV during epoch 3). Response ma- 
trices from 1998 July 9 were obtained from the 
HEASARC archive 3 . At various times 3-5 PCUs 
were active. Background estimates are derived 
from the Very Large Event (VLE) model to ac- 
count for cosmic events, internal particle genera- 
tion and South Atlantic Anomaly activation 4 . 

3 http://xte. gsfc.nasa.gov 

4 http://lheawww. gsfc.nasa.gov/ stark/pca/pcabackest.html 
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Fig. 2. — The 2-28 keV light curve of Her X-l, spanning approximately the peak of the short-high state to 
its end. Dashed lines correspond to the beginning and end of X-ray eclipse. The lower panels display the 
hardness ratios between the bands 2-5 keV, 5-8 keV and 8-20 keV. Data used to extract spin pulse profiles 
(Sec. 6) were taken from time sequences 1 to 6. 
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3. LIGHT CURVES 

Standard2 data from each visit were filtered 
to reject pointings closer than 10 degrees to the 
Earth's limb and those off-axis by 0.02 degrees or 
greater. Events were summed across pulse-height 
channels to provide a light curve in the energy 
range 2-28 keV with 100-s sampling and over all 
PCA layers and columns. This is presented in 
Fig. 2 with hardness ratios of 5-8 keV / 2-5 keV 
and 8-20 keV / 5-8 keV. Background models have 
been subtracted and data scaled by the number of 
active PCUs. 

X-ray count rates display a general trend of de- 
cline between the first and last visits, interrupted 
by dipping episodes. The event centered at MJD 
51004.7 is an eclipse of the X-ray source by the 
companion star. The same eclipse observed in the 
broad UV emission lines by HST provided con- 
straints on the location of dynamic gas in the sys- 
tem (Boroson et al. 2000). We observe residual 
counting events at mid-eclipse with a weighted- 
average rate of 1.0 ± 0.2 countss -1 PCU -1 , deter- 
mined from visits 9 and 10. From time-dependent 
spectral fits to GINGA LAC eclipse pointings, 
Choi ct al. (1994b) and Leahy (1995) argue for 
an extended, ionized, scattering region around 
Her X-l producing these residual photons, from 
where the mid-eclipse fluxes were found to be vari- 
able by a factor 7 from orbit to orbit. It was sug- 
gested that variability is related to the 35-d cy- 
cle, although no correlation has been searched for. 
However, by folding spectral models of visits 9 and 
10 (Section 4) through the response matrix of the 
LAC we find these residual counts to be consis- 
tent with the GINGA detection of 2.1 countss -1 
at <fo 5 = 0.67. 

Between MJD 51005.8-51006.8 we observe an 
extended dip with duration ~ half an orbital pe- 
riod. Although our sampling is not ideal the data 
are suggestive of dips with similar phasing and du- 
ration during both the preceding and following or- 
bital cycles. This is confirmed by the more com- 
plete sampling of the contemporaneous BeppoSAX 
observations of Oosterbroek et al. (2000). Dip in- 
gresses occur earlier each consecutive orbit, while 
the dip duration increases. Hardness ratios sug- 
gest these are, at least partially, the result of 
photoelectric absorption, although the dips re- 
main strong features at energies > lOkeV. Both 



Shakura et al. (1998) and Scott & Leahy (1999) 
folded data from the RXTE ASM archive and 
found that short-high state dips have significantly 
longer durations than main-high state dips. The 
ASM shows the width of these features is coher- 
ent from cycle to cycle over many years and there- 
fore is not related to the impending X-ray low- 
state. Long dip durations are probably the result 
of blending by numerous unresolved shorter events 
(Parmar & Reynolds 1995). 

Crosa & Boynton (1980) determined that dip 
ingresses migrate linearly towards earlier orbital 
phase over each 35-d period, where the time be- 
tween successive ingresses is 1.65 days. Although 
this is not exactly the beat between the 35-d and 
orbital cycles, Crosa and Boynton suggest that the 
dips are the result of periodic increases in the scale 
height of the outer accretion disk brought on by 
bursts of mass transfer as the L\ point of the com- 
panion star sweeps through the X-ray shadow of 
the disk. The dip period is the sum of the beat 
period and the orbital timescale of structure in 
the outer disk to travel from the accretion stream 
impact point to the pulsar line of sight. It is dif- 
ficult to see how this model incorporates dips of 
precession-coherent duration unless the accretion 
stream meets the disk at a radius that varies over 
the 35-d cycle. This would result in a range of 
orbital frequencies for the newly-arrived gas that 
broadens the dip. 

Schandl (1996) suggests the 1.65 day dip pe- 
riod results from the accretion stream skimming 
over the surface of the disk and crashing into it 
at a disk radius and azimuth that vary over the 
35-d cycle due to the disks warped shape. The 
width of the dips would most likely increase as 
the stream migrated closer to the disk's center, 
where an arbitrary surface area subtends a larger 
occulting angle for the compact object. The de- 
tailed shape of Schandl's model is not consistent 
with the inference that the stream reaches the in- 
ner accretion disk during the short-high phase of 
the cycle, but this particular disk shape also fails 
to reproduce the observed X-ray and optical vari- 
ability from the source. However, by adopting this 
type of model for the dips, it may be possible to 
constrain disk shape by measuring the duration of 
dips as a function of 35-d phase. 
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Fig. 3. — Spectral fits for each visit, where F x is the model flux, a a power-law exponent, n a the normalization 
of the power-law component, Nu a neutral absorption column density, / the partial covering fraction of the 
absorber, E Fo the energy of the Fe K feature, cr Fe the width of the iron feature, n Fo its normalization and xt 
the reduced goodness of fit measure. Error bars are 90 percent confidence limits. For these fits, / = 1 is 
fixed to represent a total covering absorber. 
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Fig. 5. — Data (crosses) and sample fits (histograms) from visits 3 (short-high state peak intensity), 9 (X-ray 
eclipse), 23 (broad dip) and 38 (short-high state tail). Residuals for each fit are in the lower panels. 



4. SPECTROSCOPY 

The current Her X-l model predicts that the 
tail of the short-high state is the result of the 
warped disk rim gradually occulting the X-ray 
source as the disk precesses on the 35-d cycle (Pct- 
terson 1977). We expect both partial-covering ab- 
sorption from the tenuous upper atmosphere of the 
disk and opaque obscuration from the thick disk at 
smaller heights above the mid-disk. Similarly, dips 
are possibly the result of periodic increases in disk 
height (Crosa & Boynton 1980), absorption from 
the "bright spot" impact site between stream and 
disk (Choi et al. 1994a), or the accretion stream 
itself with ejecta from the impact splashing above 
the disk plane (Schandl 1996). In this section we 
test whether the picture above is consistent with 
the time- varying energy spectrum of Her X-l. 

A single absorber is probably unsuitable to fit 
a system containing an accretion disk whose den- 
sity varies with vertical height, and with further 
possible absorption from winds, coronae and mag- 
netic accretion columns. However the data would 



do a poor job constraining the large number of 
free parameters a more complex fit requires. Nev- 
ertheless, we expect the partial covering model to 
provide a better statistical fit than the homoge- 
neous absorption model. We adopt a simple cut- 
off powerlaw model with cold, blanket absorption 
and compare this to a similar model with a partial- 
covering absorber. 

Under the same filtering constraints as in sec- 
tion 3, the Standard2 events were sampled with 
the full pulse-height resolution of the PCA. Chan- 
nels 80-255 (30-100 keV) were ignored as a re- 
sult of poor counting statistics. Channels 1-6 
were ignored due to uncertainties in background 
modeling. A cyclotron absorption feature of en- 
ergy 30-40 keV is usually present in the spectrum 
of Her X-l (e.g. Mihara et al. 1990; Dal Fiume 
ct al. 1998). The soft wing of this feature en- 
croaches the PCA bandpass. Including this wing 
in the model provides acceptable fits but since the 
center of the line is not sampled the line param- 
eters afford too much freedom to the low energy 
components of the model. We therefore only fit 
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data softer than 15 keV to avoid cyclotron con- 
tamination. Cutting off the energy spectrum be- 
low 3 keV avoids soft contamination from the high 
energy tails of spectral components reprocessed by 
the inner accretion disk. These were modeled as 
a 0.09 keV blackbody component and a broad Fe 
L emission complex at 0.95 keV by Dal Fiume 
ct al. (1998). 

An interstellar column density of Nn int =5.1 
xlO 19 cm -2 , consistent with the column to the 
source during the main-high state (Dal Fiume 
et al. 1998), is included in the model. The ef- 
fective absorption cross-section <r(E) was deter- 
mined assuming neutral, solar abundance material 
(Balucihska-Church & McCammon 1992). 

A broad Fe Ka emission line is detected at all 
times, except eclipse. We model this with a Gaus- 
sian of height n Fc , ct- width ct Fc and spectral en- 
ergy F Fc . Fits to simple thermal models of iron 
line plus absorbed blackbody or Bremsstrahlung 
emission proved statistically unsuitable. 

Model fitting employed xspec vlO.00. We 
consider the spectral distribution to be a power 
law, with exponent fixed at its unabsorbed value, 
a = 0.9 with normalization modeled 
by Dal Fiume et al. (1998) and approximately 
consistent with the partial-covering model used 
for the same short-high state by Oosterbroek 
et al. (2000). This is subject to a neutral absorber 
of column density Njj and partial covering frac- 
tion /. The composite model, S(E) is represented 
algebraically by: 

S(E) = e~ NH ^ E ) (f e -*H*(.E) + (1 _ f fj x 



n n E- a 



-((B-_E Fo )/2<t Fo ) 2 



V27 



(2) 



In our first model, where the absorbing column 
covers the entire X-ray source, / = 1 is fixed as 
a non- variable parameter. In a simple scenario 
where we assume the opaque structure and cool 
absorber are the lower and upper parts of a disk 
atmosphere at an arbitrary distance from the X- 
ray source, then the X-ray flux F x is expected to 
be correlated with n a , and anti-correlated with 
Nh (and / in the partial absorption models). 

Fig. 3 summarizes the best-fit parameters for 
each visit after fitting the blanket-covering model. 
The powerlaw normalization n a is correlated 



strongly with the overall flux during dips, eclipses 
and "normal" states and would indicate that much 
of the time-dependent behavior of Her X-l is the 
result of opaque obscuration rather than absorp- 
tion at soft energies. Nh is indeed anti-correlated 
with F x . Goodness of fits are illustrated in the 
lowest panel of the figure with the reduced \ 2 
statistic. Fits are adequate at the peak of the 
short-high state where there is negligible cold ab- 
sorption but their quality deteriorates during dip- 
ping events and as the high state decays. This 
verifies that a blanket-absorber is a poor model 
for explaining 35-d evolution. 

Fig. 4 summarizes the best-fit parameters for 
each visit after fitting the partial-absorber model. 
We first note that the addition of one further free 
parameter, /, has improved significantly the xt 
from fits to dip spectra, where xt ~ 1- During 
the peak of the high state, covering fractions are 
approximately with column densities consistent 
with the interstellar value. Nh is weakly cor- 
related with / and more-strongly anti-correlated 
with F x and n a . This indicates that variability in 
F x is driven by both solid-body obscuration and 
cold, partial absorption, at least within the frame- 
work of this limited spectral model. The reality is 
plausibly partial- or fully-covered absorption by a 
medium of spatially varying density, as you would 
expect from the atmosphere of an accretion disk. 

Within measurement uncertainties the energy 
of the Fe line is constant at 6.3 ±0.1 keV over 
the duration of the pointings. This is slightly 
lower than the measurement of Oosterbroek 
et al. (2000) and Her X-l observations in general. 
This could in principle result from erroneously 
adopting a = 0.9. However, adopting power- 
laws of various fixed slope provides no energy 
increase and yields either identical or poorer fit 
quality. Additionally, allowing a freedom to vary 
within the fit improves xt but the line energy 
does not increase significantly. This is identical to 
the model adopted by Oosterbroek et al. (2000), 
and there is agreement within uncertainties be- 
tween a = 0.9 ± 0.1, F x = 9xl0 10 erg s" 1 cm" 2 , 
ct Fo = 0.6 ± 0.2 keV, n Fo = (2 ± l)xl0~ 3 ph s" 1 
cm" 2 / = 1 ± 0.4 and N H = (2 ± 3)xl0 19 cm" 2 , 
for xl = 0-9. 

Line strength is correlated with the powerlaw 
component normalization. This indicates that the 
line strength is also modulated by both 35-d effects 
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and dip events. Dips have the same duration in 
both continuum and line. The inference is that 
the powerlaw component and line have a common 
locality. 

Residual counts remain during eclipse which 
can be fit with a single powerlaw of a = 1.0 ± 0.1 



and n a = (6.7 ± 2.0) x 10~ 4 ph cm" 



kcV" 1 at 



1 kcV. The addition of a line component with a 
physically suitable energy between 6-7 keV does 
not improve the fit statistically. 

Fig. 5 displays a sample of fits (the maximum 
and end of the short-high state, and both mid- 
eclipse and mid-dip). 

5. PULSE AND ORBITAL EPHEMERIDES 

A measurement of the pulse period at this 
epoch is valuable for determining the spin-up rate 
during the low-state using future data. Fig. 6 
displays the neutron star spin period history of 
Her X-l over the previous 30 years. The gen- 
eral trend is one of spin-up, but with at least 
two episodes of spin-down at <~ MJD 45000 and 
<~ MJD 49000. These spin-down episodes coin- 



cide with the two anomalous low states previ- 
ously reported by Parmar et al. (1985) and Vrtilck 
et al. (1994). 

This indicates that the anomalous low-states 
are related to a change in the location of the 
threading region, where gas is stripped from the 
accretion disk by the magnetic field of the pul- 
sar. If the characteristic distance of this region is 
closer to the neutron star than the co-rotation ra- 
dius, where the rotational velocity of the accretion 
disk exceeds the rotational velocity of the com- 
pact object, then assuming conservation of energy 
and angular momentum, the star will be spun-up. 
Conversely if the threading region occurs outside 
the co-rotation radius, we expect spin-down. Ig- 
noring any diamagnetic properties of the accret- 
ing gas (King 1993), the location of the thread- 
ing region is governed by the balance between 
magnetic pressure and the ram pressure of the 
infall. Consequently observations of spin-down 
during the anomalous low-states are naturally ex- 
plained by a decrease in the mass transfer rate 
in the inner disk. Accretion does not stop alto- 
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Fig. 7. — Times of superior conjunction of the neutron star minus the first two terms of the ^ndVrtler 01 12 23 
polynomial fit determined in Sec. 5. The dotted line is the quadratic term and the discontinuous line is two 
linear fits broken at MJD 45120. (1) The mean BATSE measurement between MJD 48100-50600. 



gether since UV and optical observations indicate 
that the companion star is still irradiated strongly 
by X-rays from the central object and inner disk 
(Vrtilek et al. 2001). The continued brightness at 
optical and UV wavelengths also seemingly rules 
out the alternative suggestion by van Kerkwijk 
et al. (1998) that spin-torque reversals are the re- 
sult of accretion disk warp flipping the inner disk 
over by 180 degrees. 

Furthermore Deeter et al. (1991) show that 
there is also an increase in the orbital period of 
the binary over time (see Fig. 7) and suggest that 
a broken linear expression provides a better fit to 
the period measurements than a continuous one. 
The break coincides with the low-state at ~ MJD 
45000, suggesting that the anomalous lows also 
have a measurable effect on the orbital rate. This 
result is less convincing than stellar spin-up mea- 
surements and requires confirmation from orbital 
timing during, and after, other anomalous low- 
states. Therefore a measurement of the orbital 
period at this epoch is equally useful. 

For each visit, GoodXenon PCA data were re- 



duced to light curves with a sample rate of 0.02- 
s over the entire detector energy response. Con- 
temporaneous background models were subtracted 
and counts scaled by the number of active PCUs. 
Data with an elevation above the Earth's limb less 
than 10 degrees or offset from the optical axis by 
greater than 0.02 degrees were ignored. Visits 7- 
11 during X-ray eclipse were also discarded. Times 
were corrected to the solar system barycenter. The 
accuracy of the on-board clock is 100 /isd -1 , al- 
though constant health checks adjust the clock fre- 
quency. Therefore stability on hourly timescales 
may be somewhat poorer. 

5.1. Orbital Pulse Delays 

In order to measure the spin period, we must 
first determine time-delays due to the motion of 
the neutron star. Data were separated into 200-s 
intervals and folded on a period close to the ex- 
pected pulse period during this epoch, -P S p m = 
1.237730-s (extrapolated by eye from the BATSE 
timings of Wilson, Scott & Finger 1997). A Gaus- 
sian function of pulse phase was fit to each folded 
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Fig. 8. — Times of pulse maximum arrival relative to the center of mass. The 2nd order orbital cphcmcris 
of Stelzer ct al. (1997) is overlayed on the data. Residuals plotted in the lower panel are accompanied by a 
linear fit. 



set and the centroid of each fit adopted as pulse 
maximum. Due to the 35-d cycle, pulse profiles 
vary over the duration of these visits and conse- 
quently this method was more practical than con- 
volving individual intervals with each other. This 
provided the fractional component of the orbital 
pulse delay. The integer component was inferred 
by direct comparison with the orbital fit of Stelzer 
et al. (1997). The orbital fit and the inferred de- 
lays are presented in Fig. 8. We find a linear trend 
over time in the residuals where the discrepancy 
between the true and assumed pulse periods re- 
sults in a gradient of -0.219 ± 0.008-sd- 1 . The 
prediction is that Fourier analysis will reveal a 
pulse period of P spin = 1.2377257 ± 0.0000004- 
s after removal of pulse delays, where the uncer- 
tainty is the 90 percent confidence limit. 

The orbital period P orb was determined by mini- 
mizing \ 2 between pulse-arrival times and a linear- 
plus-circular function (i.e. Still et al. 1994), deter- 
mining the phasing of neutron star superior con- 
junction, and pulse-arrival semi-amplitude, 
a sin i. Best fit parameters are provided in Table 1. 



5.2. Orbital Ephemeris 

was combined with the previous measure- 
ments of Deeter, Boynton & Pravdo (1981), Dccter 
ct al. (1991), Wilson, Scott & Finger (1997) and 
Stelzer et al. (1997). We fit these first with a 
2nd order polynomial to determine a new orbital 
ephemeris, provided in Table 2, where 

T ocl = T n/2 + P orb E + ^P OIb P orb E 2 (3) 

T ocI are the times of mid-eclipse and E the integer 
orbital cycle number since T v / 2 - 

The quality of fit is slightly improved if we 
adopt the suggestion by Deeter ct al. (1991) that 
the distribution is described by a broken linear ex- 
pression. The corresponding ephemerides are also 
given in Table 2 where the break is found to oc- 
cur at MJD 45120. Both fits are overlayed on the 
data in the O — C (Observed — Computed) dia- 
gram of Fig. 7, where the first two terms of the 
quadratic ephemeris have been subtracted from 
data and fits. 
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Table 1 

System parameters determined from the current pulse timings. T ff / 2 IS defined as when 

THE MEAN LONGITUDE EQUALS ir/2 (i.E. SUPERIOR CONJUNCTION OF THE NEUTRON STAR), a IS THE 
SEPARATION OF THE PULSAR AND THE BINARY CENTER OF MASS, AND i THE ORBITAL INCLINATION. 



T n/2 MJD 51004.729581(9) 
asini 13.1902(9)-s 
P orb 1.7002(3)-d 
P pin 1.237727(l)-s 



Table 2 

Best quadratic and linear fit parameters for the orbital ephemeris of Her X-l. The 

BREAK BETWEEN THE TWO LINEAR FITS OCCURS AT MJD 45120. 





T n/2 (MJD - 51000) 


P orb (d) 


P orb (dyr- 1 ) 


xl (d.o.f) 


Quadratic 


4.729549(7) 


1.700167427(9) 


-1.33(7)xl0" 8 


2.18 (19) 


Linear 1 


4.730571(50) 


1.700167790(10) 




1.13 (14) 


Linear 2 


4.729581(7) 


1.700167504(7) 




1.16 (4) 



5.3. Pulse Period 

Orbital motion was removed from the filtered 
GoodXcnon timings by subtracting a circular 
function of semi-amplitude a sin i and phased by 
Linear ephemeris 2 from Table 2. A Lomb-Scarglc 
period search was performed over the entire sam- 
ple (Press et al. 1992) and the spin period deter- 
mined as P apin = 1.237727 ± 0.000001, consistent 
with the prediction from Sec. 5.1. The error is the 
1-ct width of a Gaussian fit to the power peak in 
the Lomb-Scargle statistic distribution. 

Fig. 6 displays the appended pulse history 
where we see a general trend of spin- up, with at 
least two epochs of spin-down which coincide with 
extended X-ray low states (Vrtilek et al. 1994). 
We note that the first well-sampled episode of 
spin-down at MJD 45000 coincides with the break 
between the two linear orbital ephemerides from 
Sec. 5.2. P spin time-derivatives during the three 
epochs of spin-up were determined by linear SVD 
fits and listed in Table 3. 



6. PULSE PROFILES 

Finally in this paper we look at the evolution 
of spin pulses over the short-high phase. Both 
Deeter et al. (1998) and Scott, Leahy and Wilson 
(2000; hereafter SLW) looked in detail at the pro- 
file and spectral variability of the short-high state 
pulses. With this current epoch of data we are 
provided with short-high state sampling improved 
over those previous observations. 

We filter the GoodXcnon events into three en- 
ergy bands: 2-5 keV, 5-8 keV and 8-20 keV, 
remove pulse-delay effects due to orbital motion 
(Sec. 5.1) and extract events from six separate 
time ranges (see Fig. 2): 1: MJD 51003.8-4.2 
(peak of the short-high state), 2: MJD 51004.7- 

4.8 (mid-eclipse), 3: MJD 51004.9-5.2 (post- 
eclipse/prc-dip) 4: MJD 51006.0-6.6 (dip), 5: 
MJD 51006.8-7.0 (post-dip) and 6: MJD 51007.5- 

7.9 (the end-points of the short- high state). These 
are presented in Fig. 9. Fig. 10 presents a smooth 
greyscale representation of the pulses throughout 
the decline of the high state. In this case events 
have been filtered into 40 time bins and 30 phase 
bins. Over each time bin the mean count rate 
has been subtracted from the data. Bins which 
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Fig. 9. — Pulse profiles sampled over five time ranges (see Fig. 2). 1: MJD 51003.8 4.2 (peak of the short- 
high state), 2: MJD 51004.7-4.8 (mid-eclipse), 3: MJD 51004.9-5.2 (post-eclipse), 4: MJD 51006.0-6.6 (dip), 
5: MJD 51006.8-7.0 (post-dip) and 6: MJD 51007.5-7.9 (the end-points of the short-high state), and three 
energy bands: 2-5 keV, 5-8 keV and 8-20 keV. A and D refer to pulse components identified by Scott, 
Leahy & Wilson (2000). 
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Table 3 

Linear fits to the three episodes of spin up shown in Fig. 6. 



Epoch (MJD) 


P spin (Msy- 1 ) 


Xl (d-o.f.) 


41399-44000 


-3.81 ± 0.16 


1.39 (4) 


45761-48552 


-6.67 ± 0.08 


1.94 (22) 


49210-51006 


-4.35 ± 0.27 


2.41 (3) 



contain no data are represented by white horizon- 
tal strips. Fig. 10 shows a small phase drift in 
the pulse, modulated on the orbital period with 
amplitude a few pixels. This is most likely caused 
by a small error (within measurement uncertain- 
ties) in the orbital time delay fit. Consequently 
absolute alignment of features at different times is 
ambiguous and caution is required when aligning 
narrow features in the pulse profile. 

Using GINGA (Deeter et al. 1998) and RXTE 
observations, SLW developed a geometric model 
of the X-ray source based on the evolution of the 
pulse profile through the main- high and short-high 
states. Given an unobscured view of the accre- 
tion region, SLW argue that back-scattering off the 
two magnetic accretion curtains dominates over 
forward-scattering. However, as the cycle pro- 
ceeds, the edge of the accretion disk passes across 
our line of sight to the various emission regions, re- 
moving individual components from the pulse pro- 
file. SLW argue that emission from the far curtain 
(components B and E, in their terminology) and 
emission from the near pole (component C) are ob- 
scured by the disk at the current 35-d phases (the 
far pole is never visible in the SLW model). We 
see no obvious signatures of components B, C or E 
in the observed profiles, consistent both with this 
model and in agreement with the GINGA satel- 
lite observations of the short-high state (Deeter 
et al. 1998). 

The top three profiles of Fig. 9, extracted dur- 
ing the short-high state maximum, are consistent 
with the profiles presented by SLW from a simi- 
lar 35-d phase, showing a narrow peak at 4> Bpin = 
0.5, superimposed on a broad asymmetric modu- 
lation with minimum at spin = 0.7 and maximum 
at </> spin = 0.3, thus rising slower than it falls. Be- 
tween the D and A peaks is a minimum at Bpin = 



0.6. It is reasonable to associate the narrow peak 
with component A, which SLW model as emission 
back-scattered off the near curtain, close to the 
magnetic pole on the surface of the pulsar. We as- 
sociate the broad peak with component D, which 
is emission from the same curtain, but significantly 
higher above the surface of the star. Component 
A is harder than component D. Component A in- 
tensity and hardness ratio are directly correlated, 
while component D intensity and hardness ratio do 
not appear to be correlated over the whole pulse 
cycle. The hardness ratio has a sharp maximum 
at 4> Bpin — 0.5, corresponding to the sharp intensity 
peak A, and a broad roughly sinusoidal modula- 
tion peaking at </> spin = 0.5, half way up the inten- 
sity rise of component D. The dip at spin = 0.6 is 
not accompanied by an increase in hardness ratio, 
as expected from photoelectric absorption. 

The next row of profiles are sampled during 
mid-eclipse. During mid-eclipse the count rate is 
100 times less than the high state peak however 
hardness ratios are essentially unchanged, consis- 
tent with the findings of Choi et al. (1994). Pulse 
fractions decrease during eclipse and there is no 
compelling evidence for pulses. Consequently the 
pulses are more deeply eclipsed than the mean 
light. The residual is thought to result from a 
corona surrounding the accretion disk and neu- 
tron star and/or the companion star (e.g., Leahy 
1995). 

The third set of pulse profiles have been taken 
from the post-eclipse recovery. In the time be- 
tween the first and third pulses, component A has 
disappeared from the profiles. This is consistent 
both with the GINGA observations in the latter 
stages of the short-high state and the accretion 
disk occulting the region responsible for this com- 
ponent. Component D is slightly harder than be- 
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fore eclipse, suggesting a higher column density in 
front of this region. 

Set 4 of pulse profiles was extracted from the 
center of the next absorption dip. The hard- 
ness ratio has increased by a factor 3 over pre- 
dip pulses, expected from photoelectric absorp- 
tion, and the dip at <p spin = 0.6 has disappeared. 
Here there are two peaks at apin = 0.2 and apin 
= 0.5 and there is tentative evidence that both 
peaks are accompanied by corresponding increases 
in hardness. Unlike previous visits, intensity and 
hardness ratio are directly correlated. No pulses 
are detected in the 2-5 keV band. 

Pulses extracted after the dip and displayed as 
set 5 show a predictable decrease in hardness ra- 
tio. However the correlation between intensity and 
hardness ratio remains, as do the two peaks at 4> spia 
= 0.2 and </> apl „ = 0.5. 

In the final set of light curves, the mean in- 
tensities are a factor 10 smaller than those in the 
top row, with hardness ratios indicating a larger 
absorption column. Evidence for pulsations is 
marginal however, indicating that the pulse frac- 
tion is reduced by at least a factor 20. 

The model of SLW predicts that the pulse pro- 
file should become simpler as the short-high state 
evolves into the low state. This is the result of 
a number of emission regions being occulted by 
the disk, one-by-one as the 35-d cycle proceeds. 
This is not observed. The complexity of the pulse 
profile increases at the tail-end of the short-high 
state. Since it is not possible for component A to 
re-emerge from behind the accretion disk at some 
time during MJD 51005.2-5.6, the simplest inter- 
pretation for the maximum at </> apin = 0.5 is that 
the pulse shape of component D evolves from a 
single asymmetric peak to a double-peaked form. 
Since the hardness ratios indicate no obvious pho- 
toelectric absorption event between the two peaks, 
perhaps the easiest way to achieve this involves a 
physical change in the geometry of the accretion 
flow. This is a reasonable solution since the lo- 
cation of the threading region between disk and 
curtain will be modified cyclically with the pre- 
cession of the accretion disk. 

7. CONCLUSIONS 

We have presented 3-30 keV RXTE PCA data 
of what could arguably be one of the last normal 



short- high states of the X-ray pulsar Her X-l be- 
fore it entered an anomalous low-state of activity. 
Fig. 1 illustrates that the RXTE ASM detected 
its last major main-high state before the anoma- 
lous low at MJD 51195. There were very few (if 
perhaps no) short-high states detected after MJD 
51004, although we note Her X-l is a source under- 
going deep eclipses and dip events and the ASM 
data does not have ideal sampling. Since accre- 
tion disk structure is projected closer to the X- 
ray source during short-high states it is likely that 
the disappearance of short-high states may be a 
precursor indicating an impending anomalous low- 
state. We suggest that the event that caused the 
anomalous low may have began some time before 
it was noticed in the ASM. A study of pulse tim- 
ings over this period may reveal that spin-down 
began (or that spin-up was decelerating) long be- 
fore the high states disappeared. An interesting 
test to perform would be the measurement of P apin 
more uniformly over the epoch covered by Fig. 1 
using timings from the BATSE experiment that 
was on board CGRO (e.g., Wilson, Scott & Fin- 
ger 1997). Comparing the current epoch PCA 
data with the large sample of short- high states col- 
lected by the RXTE ASM and the limited number 
of previous pointed short-high state observations 
by GINGA, we find no evidence for this particular 
short-high state to be considered special. 

We have determined the orbital and pulsar spin 
periods at this epoch to facilitate measurements 
of P Bpin and P OIb during the anomalous low and 
during the next epoch of high-state activity. 

Broad-band timing verifies that the pre-eclipsc 
dips during the short-high state have durations 
longer than the characteristic durations of dips 
during main- high states (Scott & Leahy 1999). It 
appears likely that dip durations are related to 
the location of the gas stream impact with the ac- 
cretion disk, the trajectory of the stream and the 
shape of the disk at the impact point. Detailed 
hydrodynamical calculations may be able to place 
model constraints on these properties using dip 
observations. 

Spectrally, the short-high state and the dips 
are consistent with obscuration of a central X-ray 
source by a cloud of varying column density. The 
standard model of a warped accretion disk of finite 
scale height fits this picture well. 

Some modification to the recent pulse model of 
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SLW is required to explain the evolution of pulse 
profiles, from relatively simple to more complex, 
during the tail of the short-high state. We suggest 
(as do SLW) that geometric changes in the accre- 
tion curtains need to be considered to model the 
pulse profiles correctly. 

This paper employed RXTE All Sky Monitor 
results made publicly available by the ASM/ RXTE 
Teams at MIT and NASA/Goddard Space Flight 
Center. This work was partially funded by NASA 
grants NAG5-6711 and NAG5-7333. KOB, KH 
and HQ acknowledge research grant support from 
the UK Particle Physics and Astronomy Research 
Council. KH acknowledges support from a Beat- 
rice Tinsley Visiting Professorship at the Univer- 
sity of Texas, Austin. We thank the referee for a 
thorough contribution. 
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